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ABSTRACT 

We present the results of a five month monitoring campaign of the local active galactic nuclei (AGN) 3C120. Observations with a 
median sampling of two days were conducted with the robotic 15cm telescope VYSOS-6 located near Cerro Armazones in Chile. 
Broad band (B,V) and narrow band (NB) filters were used in order to measure fluxes of the AGN and the H/3 broad line region (BLR) 
emission line. The NB flux is constituted by about 50% continuum and 50% H/? emission line. To disentangle line and continuum 
flux, a synthetic H/? light curve was created by subtracting a scaled V-band light curve from the NB light curve. Here we show that the 
Hy6 emission line responds to continuum variations with a rest frame lag of 23.6 ± 1.69 days. We estimate a virial mass of the central 
black hole M BH = 57 ± 27 • 10 6 Af , by combining the obtained lag with the velocity dispersion of a single contemporaneous spectrum. 
Using the flux variation gradient (FVG) method, we determined the host galaxy subtracted rest frame 5 100 A luminosity at the time 
of our monitoring campaign with an uncertainty of 10% (Lagn = (6.94 ± 0.71) x 10 43 erg s~ l ). Compared with recent spectroscopic 
reverberation results, 3C120 shifts in the Rblr - Lagn diagram remarkably close to the theoretically expected relation of R oc L 05 . Our 
results demonstrate the performance of photometric AGN reverberation mapping, in particular for efficiently determining the BLR 
size and the AGN luminosity. 
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1. Introduction 



Reverberation mapping (Blandford & McKee 1982| l, where 
spectroscopic monitoring is used to measure the response delay 
t of the broad emission lines to nuclear continuum variations, 
has proven to be a powerful tool to measure the average distance 
of the BLR clouds to the central source Rblr = t ■ c. 

The spectroscopy provides us also with a velocity dispersion 
cr v of the emitting gas. Then, adopting Keplerian motion, one 
may estimate the enclosed mass, dominated by the supermassive 
black hole (Pet erson et al.|2004| and references therein). 



From theoretical considerations (Netzer 1990), the relation- 
ship Rblr k L^ gn , between H/3 BLR size and nuclear luminosity 
(5100A), should have a = 0.5. This has been investigated inten- 
sively (|Ko ratkar & Gaskell|1991[ |Kaspi et al.|2000] [2003| |Bentz| 
|et al.|2006j|2009a} , with the latest result being a = 0.519+^. 

Notably a new technique for distance determination based on 
dust-reverberation mapping has been presented (Yoshii 2002). 
Modifying these concepts, it has recently been proposed that the 
Rblr k ^agn relationship can be used as an alternative lumi- 
nosity distance indicator. The intrinsic luminosity of the AGN 
(Lagn) can be inferred from the radius of the BLR (/?blrX result- 
ing in the determination of the AGN distance. Moreover, due to 
the large luminosity and the extensive range of redshift at which 
the AGNs can be observed, the Rblr - Lagn relationship offers 
the opportunity to discriminate between different cosmologies 



and to probe the dark energy (|Haas et aL]|201 1| [Watson et al. 



201 1 1. However, these methods require that the large scatter of 
the current Rblr - Lagn relation can be reduced significantly, 



i.e. by factors up to 10, and that reverberation mapping of large 
samples can be performed efficiently. 

Recently, Haas et al. (2011 1 have revisited photometric 
reverberation mapping. They demonstrated for Arkl20 and 
PG0003+199 (Mrk335) that this method is very efficient and 
even applicable using very small telescopes. Broad filters are 
used to measure the triggering continuum, while suitable narrow 
band filters catch the emission line response. 

The estimation of the host-subtracted nuclear luminosity 
Lagn is challenging as well. One may use high-resolution imag- 
ing data and model the host galaxy profile in order to disentan- 

using HST imaging). 



gle the nuclear flux (Bentz et al. 2009a 



An alternative approach is provided by the flux variation gradi- 
ent method (FVG, by |Choloniewski| 1 98 1 1 |Winkler| 1 997) . This 
method can be easily applied to monitoring data and does not 
require high spatial resolution. 

3C120 is a nearby Seyfert 1 galaxy known to be strongly 
variable in the optical, characterized by short and long term 
variations with amplitudes of up to 2 mag on a timescale of 
10 years ( |Lyutyi|[T9"79j [Webb et aL"1|1988[ ). Sub sequent studies 
showed amplitude variations of about 0.4 mag (Winkler 1997) 
and 1 .5 mag ( |Sakata et al.|2010 1 on a timescale of 4 years mon- 
itoring. Furthermore, 3C120 was included in a complete eight 
years AGNs spectroscopic monitoring campaign conducted by 
Peterson et al. ( | 1998a >. Although 3C120 was a lower priority 



object (with 50 days average intervals between observations) the 
spectroscopic reverberation mapping results have shown that the 
H/3 emission line time response allows one to determine the BLR 



size with an uncertainty of about 50% (Peterson et al. 



Through 3C120's redshift of z = 0.0331, the H/3 line fa 



2004). 
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Table 1. Characteristics of 3C120 



a (2000)" > 


6 (2000) (1) z (1) 


Df 
(Mpc) 


H/3/ag (3) 
(days) 


(M ) 


,4 <4) 
(mag) 


Ay 

(mag) 


04:33:11.095 


+05:21:15.620 0.03301 


145.0 


38.5 


55.5 ± 26.9 x 10" 


1.283 


0.986 



(I) 



Values from NED database., (2) |Bentz et aLl|2009bf, (3) |Peterson etaT](|2004|>, (4) |Schlegel etaT1(l998). 
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Fig. 1. CAFOS spectrum of 3C120 obtained on Oct.27"' 2009. 
For illustration, the band passes of the filters used for the photo- 
metric monitoring are shaded (blue B-band, green V-band). The 
transmission of the narrow band ([OIII] 5007 A filter) is indicated 
by the red dotted line. While narrow band catches the redshifted 
Hy6 line, its flux is composed of about 50% H/3 line (red shaded) 
and 50% the continuum (black shaded). Note that for actual flux 
calculations the filter curves are convolved with the transmission 
curve of the Alta U16 CCD camera. 



the OIII filter at 5007A. Consequently, it is an ideal candidate 
for our instrumentation, a robotic 15 cm refractor. 

Here we present new measurements of the BLR size, host- 
subtracted AGN luminosity and black hole mass based on a well- 
sampled photometric reverberation mapping campaign, allowing 
us to revisit the position of 3C120 in the BLR size - luminosity 
relationship. As a lucky coincidence, Grier et al. (2012) have 
carefully monitored 3C120 spectroscopically one year after our 
campaign allowing us to directly compare our photometric mon- 
itoring results with their spectroscopic results. 

2. Observations and data reduction 

The photometric monitoring campaign was conducted between 
October 2009 and March 2010 using the robotic 15 cm VYSOS- 
6 telescope of the Universitatssternwarte Bochum, located near 
Cerro Armazones in Chilefl 

The images were reduced using IRAlQ packages and cus- 
tom written tools, following the standard procedures for image 
reduction. Light curves were obtained with a median sampling 
of 2 days in the B-band (Johnson band pass = 4330 + 500 A), 

1 More information about the telescope and the instrument has been 



published bylHaas et al. 1 201 1 



2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 




50 100 150 

Days after MJD 55105 (1.10.2009) 

Fig. 2. Observed light curves for 3C120 between Oct. 2009 and 
March 2010. The H/3 light curve is computed by subtracting a 
scaled V curve from the NB curve and re-normalizing it to mean 
= 1 . The light curves are vertically shifted by multiples of 0.2 
for clarity. The data gap between the end of December 2009 and 
the beginning of January 2010 is mostly due to strong wind pre- 
venting observations. 



V-band (Johnson band pass = 5500 ± 500 A) and the redshifted 
H/3 (NB = 5007 ± 30 A) line. Photometry was performed using 
a 775 aperture. The light curves are calculated relative to ~ 20 
nearby non-variables reference stars located on the same field, 
having similar brightness as the AGN. The absolute calibration 
was obtained using the measured fluxes of reference stars from 
SA095 field (Landolt 2009) observed on the same nights as the 
AGN. 

Additionally, we obtained an single epoch spectrum using 
the Calar Alto Faint Object Spectrograph (CAFOS) instrument 
at the 2.2 m telescope on Calar Alto observatory, Spain. The 
spectrograph's slit width was 1'.'54. The reduced spectrum is 
shown in Figure [T] 3C120 lies at redshift z = 0.0331 so that 
the H/3 line falls into the NB 5007 + 30A filter. The NB fil- 
ter effectively covers the line between velocities -3800km/s and 
+2200km/s. We determined that about 95% of the line flux is 
contained in the NB filter (red shaded line in Figure [TJ through 
line profile deconvolution. 

The characteristics of the source and a summary of the pho- 
tometry results are listed in Table [T] and Table [2] respectively. 
Note that the magnitudes and the mean of the total flux have 
been corrected for galactic foreground extinction according to 
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Fig. 3. Flux-flux diagram for the NB (OIII) and V filter. Black 
dots denote the measurement pair of each night. The red and 
green lines represent the average flux in the OIII and V band 
respectively. Fluxes were measured using circular 7.5" apertures. 
The data are as observed and not corrected for extinction. 



|Schlegel et al.| l [l998) . Observed fluxes in all bands are listed in 
Table 5. 
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3. Results and discussion 

3.1. Light curves and BLR size 

The light curves of 3C120 are shown in Figure [2] The B-band 
shows a gradual increase from the beginning of October to the 
beginning of December by 30%. Afterwards, the flux under- 
goes an abrupt drop by about 20% until mid-December 2009. 
Between the end of January and early February, the variability is 
more regular and the flux increases again to reach a third maxi- 
mum at the end of February 2010. In contrast to the steep B band 
flux increase, the narrow band (NB) flux increase is stretched 
until December 2009 to its first maximum. The sharp B band 
flux decrease in December 2009, is reflected in the NB in early 
January 2010. Thus, at a first glance, the time delay of the HB 
line against the continuum variations is 20 - 30 days. 

As usual, the precise lag determination is done via cross cor- 
relation of light curves. However, as pointed out by Haas et al. 
(201 1), the NB contains not only the emission line flux, but also 
a contribution of the varying continuum. Then cross correlat- 
ing B and NB, in principle, will result in two peaks, one peak 
from the emission line lag and one peak at zero lag from the 
auto-correlation of the continuum. These two peaks can only 
be discerned, if their separation is larger than the width of the 
auto-correlation. As illustrated in the top panel of Figure |4] in 
our case of 3C120 the two peaks are not separated, rather the 
B/NB cross correlation shows only a very broad distribution. 
As proposed by Chelouche & Daniel ( 2012| l, one could try to 
disentangle the peaks in the correlation domain, for instance by 
subtracting a B/V autocorrelation from the B/NB cross correla- 
tion. But this requires a very precise normalization of the corre- 
lations. Therefore, we prefer a more straight forward approach 



Fig. 4. Cross correlation of B and NB light curves (top) and of B 
and HB light curves (bottom). The dotted lines indicate the error 
range (±l<x) around the cross correlation. The red shaded area 
marks the range used to calculate the centroid of the lag (vertical 
red dashed line). 



by directly removing the continuum contribution from the NB 
light curve before applying the correlation techniques. 

From the spectrum we estimate that about 50% continuum 
and 50% HB emission line contributes to the narrow band fil- 
ter. The spectrum, however, was obtained with a much smaller 
aperture (slit width 1'.'54) than the photometric light curves (7'.'5 
diameter). Moreover, we have only one single epoch spectrum, 
which does not allow us to easily determine the varying contin- 
uum contribution to the NB light curve. Therefore we consider 
the NB and V band light curves. Figure [3] shows for each night 
the narrow band (OIII) and V band fluxes, calibrated to mJy. If 
the flux errors of each night are sufficiently small, then in prin- 
ciple one could subtract the V band mJy light curve from the 
NB mJy light curve to obtain the HB mJy light curve. However, 
our experience shows that this results in a rather noisy HB light 
curve. Therefore, we used an average scaling of the NB and V 
band light curves, which turned out to yield good results. In fact, 
the results, i.e. the derived lags, do not depend sensitively on the 
precise choice of the average scaling factor. 

The V band flux, on average, corresponds to ~ 56% of the 
narrow band flux (Fig. [3}. Considering that the V flux comprises 
the contributions from the continuum, the HB and the OIII lines, 
our choice of 50% continuum flux in the NB is justifiable. Thus, 
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Table 2. Photometry results of 3C120 



B 

(mag) 


V 

(mag) 


OIII 

(mag) 


f d(D 
J "total 

(mJy) 


fyW 

J tola} 

(mJy) 


foniZ 

(mJy) 


14.27-14.66 


13.92-14.18 


13.09-13.42 


7.02±0.11 


9.27±0.10 


18.03±0.26 



fBtotai , fV lo t a i and fOHI, ota i refer to the mean of the total flux ranges during our monitoring. 




4 6 
bin size (days) 



Fig. 5. Lag determined by the DCF versus bin size. The square 
symbol indicates the value used in this work (bin size of 2 days) 
and the dashed line represents the average value of 24.85. The 
data points are here connected by a line to guide the eye. 



we computed a synthetic H/3 light curve by subtracting a scaled 
V curve from the NB curve, i.e. H/3 = NB - 0.5 V. 

We used the discrete correlation function (DCF, Edels on &| 
Krohk|1988[ l to cross correlate the light curves. The cross corre- 
lation of B-band and H/3 yields a time delay of 25. 1 days defined 
by the centroid T cenl (Fig.|4|ij^] 

As usual we adopted the median sampling value (2 days) for 
the bin size in the DCF. In addition we checked whether the lag 
depends on the bin size. As discussed in detail by Rodriguez- 
Pascual et al. (1989), such dependencies could be caused by the 
scatter of the points in the light curves (bin size smaller than 
the median sampling) or by systematic noise structures (bin size 
grater than the median sampling). Figure 5] shows the lag as a 
function of bin size. Any deviation is clear y less than 2%, argu- 
ing in favour of our adopted bin size of 2 days. 

To determine the uncertainty in the time delay we applied 
the flux randomization and random subset selection method 
(FR/RSS, |Peterson et al.|1998b||Peterson et al.|2004] >. From the 
observed light curves we create 2000 randomly selected subset 
light curves, each containing 63% of the original dataj^The flux 
value of each data point was randomly altered consistent with 
its (normal-distributed) measurement error. We calculated the 
discrete correlation function for the 2000 pairs of subset light 
curves and the corresponding centroid (Fig.[6]i. This yields a me- 



3 While a correlation level of r > 0.8r max is widely used in spec- 
troscopic reverberation measurements, well sampled data allow to use 
also lower values of r (see Appendix in Peterson et al. 2004). We here 
calculated the DCF centroid above the correlation level at r > 0.5r,„ ax , 
finding that this yields more robust results for our data of 3C120, pre- 
sumably because the DCF is quite broad. 

4 Considering the Poisson probability of not selecting any particular 
point, reduces the size of the selected sample by a factor about 1/e, 
resulting in 63% of the original data (Peterson et al. 2004} . 



150 



o 100 




25 30 
Time Lag [days] 

Fig. 6. Results of the lag error analysis. The histogram shows 
the distribution of the centroid lag obtained by cross correlating 
2000 flux randomized and randomly selected subset light curves 
(FR/RSS method). The red area marks the 68% confidence range 
used to calculate the errors of the centroid. 



dian lag T ce „, — 24.4 days. Note the small lag uncertainty of 
less than 10%. Correcting for the time dilation factor we ob 
tain a rest frame lag T rest - 23.6 +1.7 days. This lag is some 
what smaller than the lag T resr 



= 27+1 days obtained by Grier 
et al. ( 2012| l, from their one year later spectroscopic reverbera 



tion mapping campaign. We will see below, that the lag differ- 
ence nicely fits to the luminosity difference of 3C120 at the two 
campaigns (Sect. 3.1.3 1. 



3.1 .1 . Virial mass of the central black hole 

Using the velocity dispersion of the emitting gas together with 
the BLR size, it is possible to determine the mass of the central 
black hole following the virial theorem: 



M BH = f- 



R-ai, 



(1) 



where o> is the emission-line velocity dispersion (assuming 
Keplerian orbits of the BLR clouds), R — c-t is the BLR size and 
the factor / depends on the geometry and kinematics of the BLR 
( |Peterson et al.|2004| and references therein). Most of the results 
presented in previous reverberation studies have been carried out 
considering only the virial product cT(Ty/G, ie. assuming a scal- 
ing factor / = 1 . 

The first empirical calibration of the factor / was obtained 
by Onken et al. (2004) who determined an average value of 



5.5+1.8, assuming that AGNs and quiescent galaxies follow the 
same Mbh-c* relationship. On the other hand, studies by |Collin] 
et al. (2006 1 show that the H/3 profile is narrower in the rms 
spectra than in the mean spectra. Therefore, a higher velocity 
dispersion is expected in the mean spectra. Hence, they suggest 
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narrow component of the H/3 profile. Under the assumption of 
standard emission-line ratios, we model this contribution (fol- 
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Fig. 7. CAFOS spectrum of 3C120, zoomed onto the H/3 line. 
The solid black curve represents the spectrum after subtracting 
the narrow [OIII]/1/14959,5007 emission lines (red dashed lines). 
The four dashed red lines (at the bottom of the H/3 profile) show 
the narrow H/3 models, for four intensities between 7% and 10% 
of the [OIII]/15007 line. The four solid blue lines represent the 
spectrum after subtracting the narrow H/3 models. Note that the 
broad faint Hell/14686 emission line can be clearly seen, but its 
possible contribution to the H/3 profile is negligible. 



a scaling factor / = 3.85 (his table 2). Subsequents studies by 
Graham e t al.| ( |20lTj ) (based on the method established by |Onken 
et al.|2004) > suggest a different scenario, reducing this first value 
to half (2.8+0.6). Therefore, the determination of the factor / 
is an open issue involving theoretical models and empirical rela- 
tionsh.png. Nevertheless, we here calculated the black hole mass 
(Mbh) adopting the value obtained by Onken et al. ( 2004 1, in or- 



der to compare our results with the literature. 

The velocity dispersion is determined from the width of the 
emission-line. It can be characterized by the second moment 
of the line profile (<Tu ne ) or by its full width at half maximum 
(FWHM). The first appears to be the most robust against mea- 
surement uncertainties ( |Peterson et al.|2 004 ) and is widely used 
in subsequent reverberation studies. Usually - for spectroscopic 
monitoring data - the line dispersion o"ii ne is measured from the 
root-mean-square (rms) spectrum, in which the contribution of 
the narrow emission lines components disappears. 

The feasibility of the use of single-epoch (SE) spectra for 
the black hole mass determination has been established in pre- 
vious investigations (e.g. [Vestergaard 20021 |Woo et al.||2007 



Penn ey et al.|2009a) . On average, uncertainties of ~30% have 
been reported for black hole mass determination from single 
epoch spectra measurements. However, the results of these stud- 
ies have detected systematic sources of uncertainties, which arise 
mainly from the emission-line velocity dispersion and narrow 



line contamination (see Denne yet al.|200 9a for details). 

We here try to determine the black hole mass of 3C120 from 
our BLR size and the contemporaneous single epoch CAFOS 
spectrum. 

The narrow emission lines [OIH]/l/i4959,5007 were mod- 
eled by single Gaussian profile and subtracted from the (SE) 
spectra (following Penn ey et al.|200 9a). We also considered the 
theoretical intensity ratio for this doublet of 2.98 determined 
by |Storey & Zeippen (|2000 ) which is consistent with obser- 
vational measurements and exhibit a negligible intensity vari- 
ation (Dimitrijevic et al. 2007). A challenge is to remove the 



lowing Penney et al.|2009a) . The [OIII]/15007 line profile was 
used to create the model of the narrow H/3 emission, adopting 
a narrow H/3 streng th between 7% a nd 10% of the [OIII]^5007 
strength (following |Woo et al. 2007). This narrow H/3 model was 
then scaled and subtracted from the observed H/3 line profile. 
Figure [7] shows the original H/3A4861 profile together with the 
narrow emission lines subtracted. 

After removing the narrow H/3 line (and the [OIII] lines), 
the H/3 profile was isolated by the subtraction of a linear contin- 
uum fit, obtained through interpolation between two continuum 
segments on either end of the line, taking into account the pos- 
sibility of local continuum contamination by Fell contribution 
and the red wing of the Hell/14686 emission line (which may 
be blended with the blue wing of H/3). Fig.|7]shows the original 
H/3/14861 profile together with the narrow emission lines sub- 
tracted. The velocity dispersion after removal of narrow lines is 
0~iine = 1504 km/ s. 

From their recent spectroscopic reverberation campaign, 
|Grier et aT7| ( |2012) > determined cr line (mean) = 1687 ±Akm/s and 
CT]me(rms) = 1514 + 65 km/ s. Note that in their mean spectra 
(their Fig. 1) the narrow H/3 component appears stronger than 
in our spectrum, suggesting that our velocity dispersion may be 
underestimated. If true, this would be one of the possible sys- 
tematic effects of single epoch spectra ( [Penney et al.|2009a) . On 
the other hand, our time lag is about 40% smaller than the result 

( 2004] ). Hence, according to photo- 



derived by Peterson et al. 



ionization physics, during our campaign the H/3 region is closer 
to the ionizing source and in consequence higher velocities of 
the gas clouds are expected. 

One shoud keep in mind that the black hole mass determina- 
tion relies on the assumption that the BLR emitting gas clouds 
are in virialized motion. Although this has been shown for a large 
sample of AGNs, there are a few exceptions that have emerged as 
a fundamental limitation for black hole mass detenriination via 
reverberation mapping (e.g. Penn ey et al.|20 09b Gaskell 2010| ). 
Single epoch spectra can not identify such different velocity sig- 
natures of H/3 clouds with a complex kinematical behavior. 

Using the previous value (with 25% uncertainty adopted and 
our rest frame time lag t = 23. 6d, the virial black hole mass is 
M V iriai = (10±5)x 1O 6 M which is consistent, inside t he margins 



Grier et al. 



of errors, with M viriaI = (12 + 1) x 10 M Q derived by 
(2012 ) from their recent spectroscopic monitoring campaign and 
with M V i ria i = (10 + 5) x 10 6 M G derived via spectroscopic re- 
verberation mapping by Peterson et al. (2004). Considering the 
factor / = 5.5 ± 1.8, we determine a central black hole mass 
M BH = (57 + 27) x 1O 6 M . 



3.1.2. Host-subtracted AGN luminosity 

In order to determine the pure AGN luminosity, commonly at 
5 100 A, the contribution of the host galaxy to the nuclear flux 
has to be subtracted. The contribution of the host galaxy to the 
nuclear flux of 3C120 has been studied by Bentz et al. (2006, 
2009a) using high resolution HST imaging as well as by|Winkler| 



eTaLl ( [T992] ), |Winkler| ^997} and |Sakata et al!) ( |2010| ) using 



the so-called flux variation gradient method, which was origi- 
nally proposed by Choloniewski ( 1981 ). We here apply the FVG 
method to our 3C120 data and compare our results with the pre- 
vious ones. Because the FVG method appears to be not widely 
known, we start with some comprehensive explanations here. 
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3.1.2.1 On the flux variation gradient method 

For multi-band AGN monitoring data, Choloniewski noticed 
a linear relationship between the fluxes obtained in two differ- 
ents filters taken at different epochs. Using this relationship he 
demonstrated the existence of two component structure to the 
spectra of Seyfert galaxies. The two component structure as- 
sumes that one component (the host) is constant in time and 
that the second component (the AGN) has a strong variability. 
Despite of this strong variability the spectral energy distribution 
(SED) of the AGN does not change. Choloniewski plotted the 
observed UVB fluxes for 40 Seyfert galaxies as 3 dimensional 
vectors. This gave a clear geometrical interpretation of the to- 
tal flux variations and allowed for the decomposition of the total 
vector flux 0(f) into a constant vector F and a variable vector 
/(f) (his figure 1). 

Following the geometrical representation from 
|Choloniewski ( 1981| l, is possible to write the two compo- 
nent structure using two different filters (/, j) as: 



10 



<f> i ,j(v,t) = F u (v)+f iJ (v,t) 



(2) 



where the observed flux Fy(v) and fij(y, t) are the constant and 
variable components in two arbitrary filters respectively. 

With the - at optical wavelengths observationally corrobo- 
rated - assumption that the shape of spectrum of the variable 
component does not change, one deduces: 



(3) 



fM _ 
mt) ~° Ji 

where ap is a constant. Using equations (2) and (3) one obtains 
the equation of a straight line in the plane (0,, 

4>j(v, t) = ajrfiiv, t) + bji (4) 

where by = F/v) - a i7 F,(v) 

Thus the coefficient ay, represents the color index of the vari- 
able component. This coefficient of proportionality has also been 
called the flux variation gradient (FVG). It is denoted by the 
symbol Y by I Winkler et al. ( 1992 1. Both, ay and bp coefficients 
have to be determined by a linear regression analysis. 

The choice of which regression method to use depends on 
the degree of knowledge about the data, especially on how well 
separated the dependent and independent variables are, the mea- 
surement errors, the intrinsic dispersion of data around the best 
fit and others factors (see |Isobe et al.|199 0, for details). 

( 1992 1 proposed a new method to separate 



Winkler et al. 



the nuclear flux from the host galaxy, based on the flux-flux di- 
agrams by |Choloniewski| ( |1981[ ). Winkler et al. monitored 35 
southern Seyfert galaxies using UBVRI multi-aperture photom- 
etry to estimate the colors of the host galaxy. They plotted the to- 
tal fluxes obtained through 20" and 30" apertures, together with 
a line representing the fluxes from a component having the same 
colors as the host galaxy. The intersection of this line with the 
linear regression fit obtained from the total flux represents then 
the contribution from the host galaxy. As the observed source 
varies in luminosity, the fluxes in the FVG diagram will follow 
a slope representing the AGN color while the host will show no 
variation. The nuclear flux is then calculated by subtracting the 
constant host galaxy component (obtained by the FVGs method) 
from the total flux. 

The studies by | Cho loniewski (19811 and Winkler et al. 
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Fig. 8. Flux variation gradient diagram of 3C120. The solid lines 
delineate the bisector regression model yielding the range of the 
AGN slope. The dashed lines indicate the range of host slopes 
determined by |Sakata et al.| ( |2010| for 1 1 nearby AGN. The in- 
tersection between the host galaxy and AGN slope (red area) 
gives the host galaxy flux in both bands. The two vertical dot- 
ted and dashed lines show the host flux determined by Bentz 
|et al.| ( [2006| and Bentz et al. (2009a) respectively. The solid line 
shows the host flux obtained by Sakat a~et al.| ( |2010"| . The dash- 
dotted blue lines represent the range of the AGN flux in both 
filters. Fluxes (black dots) were measured in a 7.5" aperture pho- 
tometry and corrected for galactic foreground extinction. 



are different but stay constant with time. Although these results 



have been corroborated by numerous authors (e.g., |Yoshii et al. 


2003; Suganuma et al. 


|2006 Doroshenko et al. 2008 Sakata 


et al.|2010|l, they are sti 
sive discussion we refer 


1 a matter of a debate. For a more exten- 
the reader to Sakata et al. 20101 



A valuable study on the determination of the nuclear flux 
contribution was conducted by Sakata et al. ( 2010| l, using both 



HST imaging and the FVG method. In this study the flux of the 
host galaxy was estimated for 11 nearby Seyfert galaxies and 
QSOs. These fluxes were measured in B, V and / bands by us- 
ing surface brightness fitting to the high resolution Hubble Space 
Telescope (HST) images and from the MAGNUM observations. 
The authors find a well defined range for the host galaxy slope 
of 0.4 < r'gy' < 0.53, which corresponds to host galaxy colors 
of 0.8 < B - V < 1.1. The results - in particular on 3C120 - 



obtained by Sakata et al. (20101 are consistent (within measure- 
ment errors) with those obtained by Winkler ( 1 997), [Bentz et al.| 
( |2006l > and |Bentz et aT1 ( |2009a| i. 

3.1.2.2 Application to 3C120 

Figure [8] shows the B and V fluxes obtained during the same 
nights and through the same aperture in a flux-flux diagram. The 
host color range is taken from Sakata et al. (2010 1 and drawn 



( 1992) observationally established at optical wavelengths like in 
the BVRI bands that the AGN and host colors, i.e. flux ratios, 



from the origin of ordinates (dashed lines). For the varying AGN 
flux (black dots), we use linear regression in order to determine 
a range of possible AGN colors (solid thin lines). The cross- 
section of the host and AGN slope allows us now to split the su- 
perposition of fluxes in both filters. Note that the B and V fluxes 
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Table 3. Linear regression results. 



Method 


a" 


a-„ 


b b 


o-b 


OLS(Y/X) 


1.04 


0.04 


-2.60 


0.53 


OLS(X/Y) 


1.21 


0.05 


-4.22 


0.33 


OLS Bisector 


1.12 


0.04 


-3.38 


0.39 


Orthogonal Regression 


1.13 


0.05 


-3.48 


0.41 


RMA 


1.22 


0.06 


-4.18 


0.49 



{a) Slope. 

<A) Intercept coefficient. 



have been corrected for galactic foreground extinction according 
to |Schlegel et"aT1 ( [T998) >. 

In both filters, the total flux (AGN+Host) contains a con- 
tribution from the emission lines originating in the narrow line 
region (NLR). However, this contribution is less than 10% of the 
host galaxy flux in the B and V band ( jSakata et aT71|2QlQ) >. We 
here define the host galaxy to include the NLR line contribution. 

Flux Variations Gradients (FVGs) were evaluated by fitting 
a straight line to the data using five methods of linear regression; 
OLS(Y/X), OLS(X/Y), OLS bisector, Orthogonal regression and 
Reduced Major-Axis (RMA) were used, depending on the cor- 
responding error treatment of each method. Our regression al- 
gorithms are based on the formulas of Isobe et al. ( 1990) (their 
Table 1) and the variance for each method was calculated using 
the same formulation. The statistics of each linear regression fit 
are listed in Table [3] 

The bisector linear regression line and the OLS(X/Y) method 
yields a linear gradient of Fbv = 1.12 ± 0.04 and Fbv = 
1.04 + 0.04 respectively. The results are consistent, within the 
uncertainties, with Ybv = 1 - 1 1 ± 0.02 determined by Sakata et al. 



p0T0| > and T BV = 1 .02 ± 0.07 obtained by |Wmkler| ( |l 997| l. While 
the uncertainties to the OLS bisector and OLS(Y/X) regression 
slopes are lower (about 4%). According to Isobe et al. (1990), 
the OLS bisector method is the most suitable in order to de- 
termine the underlying relationship between the variables and 
it was also considered in previous studies of FVG method by 
Winkler ( 1997). Therefore, we adopted the OLS bisector linear 
regression to define the range of AGN slope. 

Averaging over the intersection area between the AGN and 
the host galaxy slopes, we obtain a mean host galaxy flux of 
(2.17 ± 0.33) mJy in B and (4.58 ± 0.40) mJy in V. Our host 
galaxy flux derived with the FVG method is consistent with the 
values f B ~ 2.10 mJy and fy w 4.73 mJy obtained by Sakata 
et al.|p010) ). Note that one has to coadd the values of Tables 5 
and 8 of Sakata et al., in order to include the flux contribution of 
the narrow lines ([OIII]^4959, 5007, HB, Hy) to each filter. 

We have also compared our results with those obtained by 
Bentz et al] ( |2006) > and |Bentz et al] ( |2009a|) through mod- 
eling of the host galaxy profile (GALFIT7 |Peng et al.||2002l ) 
from the high-resolution HST images. On the nucleus-free im- 
age and through an aperture of 5" x 7'.'6, Bentz et al. (2006) 
determined a rest- frame 5100A host-flux f 510 oA =1-82 X 



10 l5 ergs l cm 2 A 1 



Using the color term factor from the HST- 
F550M filter to 5100A (their Table 3) we deduced the flux for the 
F550M filter to be F F550M =1.74 x 10~ 15 ergs' 1 cm- 2 A" 1 . This 
corresponds after exinction correction to 4.62mJy. Moreover, 
considering the contribution of the narrow lines in the V- 
band (l.235mJy), determined by Sakata et al. ( 2010| ), the pre- 
vious value translates to 5.86m7y. In a subsequent investiga- 
tion, 

~~ '-1 



ergs 



-2) 



Bentz et al.| ( |2009a[ ) determined a host galaxy flux of 
0.78 x 10~ 15 - 
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Fig. 9. Interpolated rest frame flux at 5 100 A from fB ACN and 
/Vacn fluxes (black dots) for an AGN continuum with power 
law spectral shape (F v oc v") (solid lines). 



and adding the contribution of the narrow lines yield a value 
of 2>3\mJy. The difference between the two results by Bentz 
et al. lies mainly in the type of model used for the decomposi- 
tion of the galaxy. The first study considered a general Sersic 



function for modeling bulges ( |Bentz et al. 2006). The second 
study performed the modeling with variations and improve- 
ments to the original profile (Bentz et al. 2009a i. Our value 
(fVhost = 4.58m/y) falls exactly in the middle between both 
values, simply suggesting that our determination is consistent, 
within the error margins, with those determined by Bentz et allj 

The AGN fluxes at the time of our monitoring can be deter- 
mined by subtracting the host galaxy fluxes from the total fluxes. 
During our monitoring campaign, the total B fluxes lie in the 
range between 5.81-8.39 mJy with a mean of (7.02 + 0.11) mJy. 
The total V fluxes lie in the range between 8.17-10.38 mJy with 
a mean (9.27 + 0.11) mJy. The host galaxy subtracted average 
AGN fluxes and the host galaxy flux contribution of 3C120 are 
listed in Table g] 

Also listed in Tabled are the interpolated rest frame 5 100 A 
fluxes and the monochromatic AGN luminosity /IL^agn) at 
5 100 A . The rest frame flux at 5 100 A was interpolated from 
the host-subtracted AGN fluxes in both bands, adopting for the 
interpolation that the AGN has a power law SED {F v oc v") with 
an spectral index a = \og(fB ACN /fV A GN)/ log(v fi /vy), where v B 
and vy are the effective frequencies in the B and V bands, re- 
spectively. The error was determined by interpolation between 
the ranges of the AGN fluxes +cr in both filters, respectively, as 
illustrated in Fig. [9] 

To determine the luminosity, we used the distance of 145 
Mpc ( [Bentz et al.||2009a| ). This yields L A gn = (6.94 ± 0.71) x 



which after extinction correction 



5 The actual debate on the host flux discrepancies cannot be solved 
here. We just note that 3C120 (together with another three objects) 
shows larger residuals in the surface brightness fit of the MAGNUM V- 
band images (Fig. 17 of Sakata et al. 20T0) compared to the HST images 
(Fig. 14 of |Sakata et al.|2010| and Fig. 3 of |Bentz et al.|2009a>. This may 
suggest a possible host overestimation. On the other hand, [Bentz et al. 
( 2009a I mentioned a possible host underestimation caused by the uncer- 
tain sky background. Futhermore, our aperture area (7'.' 5 in diameter) is 
only 16% larger than that of Peterson's campaign (5" x 7'.'6). This to- 
gether with an additional bandpass conversion factor between the HST- 
F550M filter and our Johnson-V filter would increase the Ben tz et al.l 
(2009a ) host flux by about 10%. This suggests that the most contempo- 
raneus value quoted by Bentz et al. ( 2009a I could be underestimated. 
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10 43 ergs . The AGN luminosity during our campaign is about 



half the mean value of (12.40 + 2.6)x I0 43 ergs 1 found by|Bentz| 



tainty in our measurements (7%) is presumably due to the 
well sampled photometric reverberation data, 
et al. ( 2009a| ) for their campaign about 15 years ago and slightly 2. The black hole mass Mbh — 57 + 27 x 1O 6 M is c onsistent 
lower than the value (9 + 1) x 10 43 ergs~' found by |Grier et al.| with the value of 55.5 + 26.9 x 1O 6 M derived by Peterson 



(2012 ) for their campaign one year after ours 



3.1.3. The BLR size - luminosity relationship 

From spectroscopic reverberation mapping, the relationship be- 
tween the H/3 BLR size and the lumi nosity (5 100 A) of the AGN 
Rm r oc L a h as been established by Kaspi et al. ( 2000 1. Bentz 
|et al.| ( [2009a| determined an improved slope of a — 0.519"'^^ 
trom several spectroscopic reverberation mapping campaigns. 
Figure [T0| shows the /?blr and Lacn values obtained by 



Bentz et al . (2009a), the result for 3C120 obtained by |Grieretal 
( |2012| ) and the two objects Arkl20 and PG0003 st udied with 
well sampled photometric reverberation mapping by Haas et al. 
( 201 l| l. While the relationship appears to be well defined, many 
objects have large uncertainties yet and/or lie quite off the regres- 
sion line. Part of this may simply be due to the fact that AGN are 
complex objects, but part of the dispersion may be due to poor 
early reverberation data. We just note that the new position of 
both Arkl20 and PG0003 lies closer to the regression line in 
Fig _ 



10 



For 3C120 we have now three positions in the R-L dia- 
gram: one from Peterson's spectroscopic monitoring campaign 
in 1998, quite on the regression line but having large uncer- 
tainties presumably due to sparse sampling; one from the spec- 
troscopic monitoring campaign in 2010/2011 by Grier et al. 
2012); and one from our photometric monitoring campaign in 



2009/2010. The last two campaigns both have a good time sam- 
pling and small uncertainties. The striking result from these two 
campaigns is that the slope between these positions of 3C120 is 
a = 0.504, hence dueing the brightness changes 3C120 moves 
parallel to the theoretically expected slope (red arrow in Fig. 10 1. 



With increasing luminosity the BLR size grows proportional to 
the square root of the luminosity, and these changes appear to oc- 
cur rather quickly, i.e. within days or weeks. However, this does 
not mean that the BLR gas clouds move that rapidly in- or out- 
wards. Baldwin et al. ( 1995 ) have presented a model of locally 
optimally emitting clouds. 

It successfully explains the emission line characteristics of 
quasar spectra. This model provides also a nice explanation for 
the R-L relation. In the BLR the gas clouds actually populate 
a large range of radii, but, depending on the density and ioni- 
sation parameter, only at a suited narrow distance range (from 
the illuminating central power house) the gas clouds efficiently 
convert the nuclear continuum radiation into line emission. In 
this picture, a central continuum brightness variation can quickly 
change the BLR size proportional to L. 



4. Summary and conclusions 

Using a robotic 15 cm telescope located at an excellent site, 
we have performed a five months monitoring campaign for the 
Seyfert 1 galaxy 3C120. We determined the broad line region 
size, the virial black hole mass and the host-subtracted AGN lu- 
minosity. The results are: 

1. The time lag j res , — 23.6 +1.69 days, obtained from cross 
correlation of the H/3 emission line with the optical contin- 
uum light curve, has changed over one decade, but the phys- 

1/2 

ical relation Rblr 00 ^agn-isioo st ^ holds. The small uncer- 



et al. ( 2004[ ) from spectroscopic reverberation mapping data 
and with the value of Mbh = 67 ± 6 x 10 6 determined by 
|Grier et aL] ( |2012| ). However, the high uncertainties in the 
black hole mass, with respect to the most contemporaneus 
result, are attributed to the adopted 25% uncertainty for the 
line velocity dispersion and also considering the uncertainty 
introduced by the scale factor /. 

Using the flux variation gradient method (FVG) and a con- 
servatively limited host galaxy color range, it is possible to 
find the host galaxy subtracted AGN luminosity of 3C120 at 
the time of our monitoring campaign to be Lagn = (6.94 ± 

'3 

ergs 



0.71) X 10 43 --— 



The new results obtained for the BLR size and AGN lumi- 
nosity of 3C120 fit well into the BLR size-luminosity diagram. 
We conclude that photometric reverberation mapping is an at- 
tractive method with the advantage to efficiently measure the 
BLR and the host-subtracted luminosities for large samples of 
quasars and AGNs. Not only applicable with small telescopes, 
photometric AGN reverberation mapping could become a key 
tool for the upcoming large monitoring campaigns, for instance 
with the Large Synoptic Survey Telescope (LSST). This could 
be an important step forward in order to constrain cosmologial 



1/2 

AGN,^5100 



relationship in order to 



parameters from the Rblr k L 
determine quasar distances and to probe the dark energy (Haas 
et al.|201 l|Watson et al.|20TT i as well as to enlarge the current 
statistics for black hole masses. 
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Table 4. Host galaxy subtracted AGN fluxes of 3C120 



fBl, ost 


f B AGN 


fVhost 


M GN 


/ AGA r((l+z)5100A 


AL^GtfSlOOA 
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(mJy) 


(mJy) 


(mJy) 


(mJy) 


(10 43 erg^') 
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4.58±0.40 


4.68±0.41 


4.72±0.40 


6.94±0.71 



(1) AGN fluxes values fB AGN = fB total - fB lwst and fV AGN = fV tota , - fV hosl , with uncertainty range cr AGN = (<rj olal + crf w J° 
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Table 5. B, V, 0111 and HyS Fluxes in mjy. 
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Table 5. continued. 
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.078 


1 


.944 


±0.151 


3.448 


± 0.032 


5 


.960 ± 


0.092 


4 


,236 ± 


0. 


093 


55234 


.121 






± - 




± - 


5. 


.928 ± 


0.094 


4 


.204 ± 


0. 


101 


55236 


.070 


1 


.966 


±0.018 


3.484 


± 0.034 


5 


.957 ± 


0.103 


4 


.215 ± 


0. 


120 


55238. 


.070 


2 


.019 


± 0.032 


3.494 


± 0.059 


6 


.027 ± 


0.089 


4 


.280 ± 


0. 


090 


55240. 


.070 


2 


.061 


± 0.028 


3.544 


±0.018 


6 


.136 ± 


0.101 


4 


.364 ± 


0. 


101 


55241. 


,051 


2 


.105 


± 0.044 


3.585 


±0.017 


6 


,174 ± 


0.106 


4 


,381 ± 


0. 


106 


55242. 


051 


2 


.114 


±0.021 


3.622 


±0.038 


6 


.281 ± 


0.102 


4 


.470 ± 


0. 


102 


55243. 


.059 


2 


.131 


±0.018 


3.635 


± 0.061 


6 


.261 ± 


0.099 


4 


.444 ± 


0. 


099 


55244. 


.039 


2 


.168 


±0.011 


3.654 


±0.015 


6 


.296 ± 


0.116 


4 


.469 ± 


0, 


116 


55245. 


.059 


2 


.159 


± 0.020 


3.643 


±0.014 


6 


.355 ± 


0.099 


4 


.533 ± 


0. 


099 


55246. 


.059 






± - 




± - 


6 


.394 ± 


0.092 


4 


.572 ± 


0. 


092 


55250. 


.051 






+ - 




+ - 


6 


.343 ± 


0.121 


4 


.523 ± 


0. 


121 


55251 


.051 


2 


.332 


±0.013 




± - 


6 


.282 ± 


0.137 


4 


.462 ± 


0. 


137 


55252. 


.051 






± - 


3.640 


±0.017 


6 


.273 ± 


0.123 


4 


.453 ± 


0. 


123 


55254. 


039 






± - 


3.608 


±0.013 


6 


,203 ± 


0.122 


4 


.399 ± 


0, 


122 


55258. 


.051 






± - 




± - 


6 


,130 ± 


0.121 


4 


,326 ± 


0. 


121 


55261. 


.051 






± - 




± - 


6 


.148 ± 


0.094 


4 


.343 ± 


0. 


094 


55262. 


039 






± - 




± - 


6 


.067 ± 


0.097 


4 


.263 ± 


0. 


097 


55263. 


.039 






± - 




± - 


6 


,232 ± 


0.063 


4 


.428 ± 


0, 


063 



